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ABSTRACT 



Context. Measurements of intracluster gas temperatures out to large radii, where much of the galaxy cluster mass resides, are impor- 
tant for using clusters for precision cosmology and for studies of cluster physics. Previous attempts to measure robust temperatures at 
cluster virial radii have failed. 

Aims. The goal of this work is to measure the temperature profile of the very relaxed symmetric galaxy cluster Abell 2204 out to large 
radii, possibly reaching the virial radius. 

Methods. Taking advantage of its low particle background due to its low-Earth orbit, Suzaku data are used to measure the outer 
temperature profile of Abell 2204. These data are combined with Chandra and XMM-Newton data of the same cluster to make the 
connection to the inner regions, unresolved by Suzaku, and to determine the smearing due to Suzaku 's point spread function. 
Results. The temperature profile of Abell 2204 is determined from ~10 kpc to ~1800 kpc, close to an estimate of rioo (the approxima- 
tion to the virial radius). The temperature rises steeply from below 4 keV in the very center up to more than 8 keV in the intermediate 
range and then decreases again to about 4 keV at the largest radii. Varying the measured particle background normalization artificially 
by ±10% does not change the results significantly. Several additional systematic effects are quantified, e.g., those due to the point 
spread function and astrophysical fore- and backgrounds. Predictions for outer temperature profiles based on hydrodynamic simula- 
tions show good agreement. In particular, we find the observed temperature profile to be slightly steeper but consistent with a drop of 
a factor of 0.6 from 0.3 r 2 oo to r 2 oo, as predicted by simulations. 

Conclusions. Intracluster gas temperature measurements up to 7200 seem feasible with Suzaku, after a careful analysis of the different 
background components and the effects of the point spread function. Such measurements now need to be performed for a statistical 
sample of clusters. The result obtained here indicates that numerical simulations capture the intracluster gas physics well in cluster 
outskirts. 

Key words. X-rays: galaxies: clusters - galaxies: clusters: individuals: Abell 2204 



1. Introduction 

Cosmologically, the most important parameter of galaxy clus- 
ters is their total gravitational mass. X-rays offer an attractive 
way to determine this mass through measurements of the 
intracluster gas temperature and density structures. X-rays are 
also a unique tool to study the physics of the hot cluster gas; 
e.g., gas temperature profiles allow constraints on (the suppres- 
sion of) heat conduction. Consequently, constraining cluster 
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Unfortunately, making these measurements is quite chal- 
lenging in outer cluster regions. Even with XMM-Newton and 
Chandra, it is very difficult to determine temperature profiles re- 
liably out to more than about 1/2 the cluster virial rad ius, r vu (see 
the references above but also lSolovveva et aT1l2007l for a differ- 
ent view). As a result, only about 1/8 of the cluster volume is 
actually probed. The primary reason is not an insufficient collect- 
ing area or spectral resolution of current instruments: the limit- 
ing factor is the high particle background. Here, the X-ray CCDs 
onboard Suzaku come into play. Owing to its low-Earth orbit and 
short focal length, the background is much lower and more sta - 
ble than for Chandra and XMM-Newton dMitsuda et al.1 12007). 
making Suzaku a very promising instrument for finally settling 
the cluster temperature profile debate. And indeed, some of us 
have recently succeeded in measuring the temperature and metal 
abundance in the outskirts of the merging clusters A399/A401 
(Fujita et al. 2008j). Here, we go one step further and also con- 
firm this great prospect for outer cluster regions whose emissiv- 
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ity is not enhanced by merging activity and determine the tem- 
perature profile of the regular cluster Abell 2204 out to ~1800 
kpc with Suzaku. This radius is close to an estimate of r2oo; i-e-, 
the radius within which the mean total density equals 200 times 
the critical density, often used as approximation to the virial ra- 
dius. Note that shortly after this paper was submitted, another 
study of a cluster temperatur e profile towards ve ry large radii 
with Suzaku was submitted bv lGeoree et al.l (120081) . 

Throughout, we assume Q m = 0.3, Qa = 0.7, and Hq = 71 
km/s/Mpc; i.e., at the redshift of Abell 2204 (z = 0.1523), 1' = 
157 kpc. 



2. Observations, Reduction, Analysis 

2.1. Chandra and XMM-Newton data 

Abell 2204 was analyzed with Chan dra as part of the 
HIFLUGCS dReiprich & Bohringerl20"02l) follow- up. The reduc- 
tion an d analysis are similar to those outlined in iHudson et alJ 
(2006). The total clean exposure of the two archival observation s 
is 18.6 ks. Further details are described in Hudson et al. (2009). 

The XMM-New ton analysis of Abell 2204 is published 
(Zha ng et al.1 l2008h . Extrapolating the mass profile using the 
published gas density and temperature profile, we found r2oo = 
11.75' = 1840 kpc. 

2.2. Suzaku data 

Abell 2204 was observed with Suzaku on September 17-18, 
2006. We started from the cleaned event files (processing version 
1.2.2.3) and reran the ftools xisputpixelquality and xispi, using 
HEAsoft 6.2 and CALDB 070409. Using xselect, the event files 
were further filtered by applying the criteria STATUS=0:65535, 
COR>6, and ELV>10. The good exposure amounts to about 
50 ks for each of the four XIS using both 3x3 and 5x5 edit- 
ing modes. Figure Q] (left) shows the combined image and the 
regions selected for spectral analysis. 

The widths of the annuli were chosen to be about twice as 
large as the half power diameter of the point spread function 
(PSF). Combined spectra were extracted using the filtered 3x3 
and 5x5 editing mode events files. Using four regions from each 
of the four detectors, this resulted in 16 spectra total from the 
source observation. The spectra were grouped to have at least 50 
counts in each energy channel. 

We checked the spectra for contamination by solar wind 
charge exchange emission. The ACE SWEPAM proton flux was 
less than 6xl0 8 protons/s/cm 2 during the observation . Therefore, 
no strong flaring is expected dFuiimoto et al.l 120071) . Moreover, 
we checked the soft band (<2 keV) XIS-1 lightcurve of the local 
background region and did not find any indication of flaring. 

Response files were created with xisrmfgen and xissimarf- 
gen. We generated the first Ancillary Response File (ARF) by 
feeding an image, constructed from the double /3 model fit to 
the Chandra surface brightness profile, into xissimarfgen. This 
ARF was later used to model the cluster emission. For the sec- 
ond ARF, which was used when modeling the Galactic fore- and 
cosmic X-ray background, we assumed a uniform photon dis- 
tribution. We simulated 10 7 photons per detector using medium 
sampling. 

Night Earth data were used to create spectra of the parti- 
cle background for each region and each detector, weighting 



the spectra by cutoff rigidity^] These particle spectra were sup- 
plied as background spectra to the XSPEC (version 12.3.1) fit- 
ting routine. The cluster emission was modeled with an ab- 
sorbed thermal model (phabs*apec). The cosmic X-ray back- 
ground and Galactic emission were accounted for with an ab- 
sorbed powerlaw and an unabsorbed thermal model, respectively 
(phabs *pow + apec) . 

All 16 spectra were fitted simultaneously. We froze the 
hydro gen column density at 6.07 x 10 20 cm" 2 (Kalberl a et al. l 
120051) . the cluster redshift at z = 0.1523 dStruble & Rood! 
1987), the redshift and metal abundance of the Galactic com- 
ponent at z — and A = 1 solar, and the powerlaw pho- 
ton index at 1.41 ( Kus hino et alj |2002). The remaining param- 
eters of the fore- and background components were left free 
and linked acr oss all regions and detec tors. The temperatures, 
metallicities dAnders & Grevessei ri989). and normalizations of 
the cluster emission were left free and linked across the four 
detectors. Annulus 4 served to help constrain the fore- and 
background components because no significant "contaminating" 
cluster emission is expected beyond the estimated r2oo> so the 
normalization of the thermal cluster model was frozen at in 
this annulus. 

In the spectral fitting for this hot cluster, we ignored all 
photons with energies <0.7 keV. We did this to minimize any 
systematic uncertainties due to the correction for the effective 
area degradation by the contamination, due to the Galactic fore- 
ground subtraction, and due to the influence of possible low in- 
tensity solar wind charge exchange emission. The upper energy 
cut used for the fitting was set at 10 keV (8 keV for XIS-1, due 
to remaining calibration uncertainties). The Si K edge (1.8-1.9 
keV) and, for XIS-0, the calibration Mn K line (5.7-6.0 keV) 
were excluded. The range 7.2-8.0 keV was ignored because 
of the presence of the strong instrumental Ni K line. Overall, 
the model provided a good description of the data, resulting in 
xl ei = 1 .20 for 4097 degrees of freedom. 

3. Results 

The spectra and best fit models are shown in Figure[2] In the left 
panel, we show all 16 spectra (data points with errors) and cor- 
responding best fit models (solid lines). The particle background 
was subtracted from the data points. Each spectrum has three 
model lines, one thin line representing the cluster emission, an- 
other thin line representing the combined Galactic fore- and cos- 
mic X-ray background, and the thick line the sum of these (i.e., 
the thick lines should match the data points). The main purpose 
of this plot is to show that there are no strong residuals anywhere. 

In the right panel, we show the same plot for the most im- 
portant region (bin 3). It is instructive to compare the role of the 
different fore- and background components relative to the clus- 
ter emission. Therefore, we show here also the particle back- 
ground spectra as stars with thick error bars (which were sub- 
tracted from the data) as well as the Galactic fore- and cosmic 
X-ray background components individually. For clarity, let's fo- 
cus on the XIS1 spectra; i.e., the red data points and lines. The 
thick line matches the data points pretty well, as it should. Now 
let's start at the low energies: The uppermost thin line repre- 
sents the Galactic foreground component. This component dom- 
inates over all other components up to about 1 keV. The next 
line represents the cluster component, which starts to dominate 

1 See http://www.astro.isas.jaxa.jp/suzaku/analysis/xis/nte/ for de- 
tails. Additionally, we discarded a very small fraction of events by re- 
quiring T_SAA_HXD>436, resulting in 795 ks total exposure time. 
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Fig. 1. Le/f panel: Suzaku image of A2204 combining all four XIS detectors. Also indicated are the regions used for the spectral 
analysis (green) and the estimated r-2100 (white). The outermost annulus (annulus 4, beyond r2oo) was used to aid in the local fore- 
and background estimation. The elliptical region contains emission from unrelated sources and was excluded from the spectral 
analysis. Right panel: Gas temperatures as function of radius as measured with Suzaku (filled circles), XMM-Newton (diamonds), 
and Chandra (triangles). For each Suzaku bin three temperatures are shown. The upper and lower temperature values were obtained 
by artificially decreasing and increasing the particle background normalization by 10%, respectively. The bins are slightly offset in 
this plot for clarity. All errors are given at the 90% confidence level. 
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Fig. 2. Left panel: Particle background subtracted Suzaku spectra from all four regions (innermost region at the top, outermost 
region at the bottom) and from all four detectors (e.g., for the innermost region: black: XIS0, red: XIS1, green: XIS2, blue: XIS3). 
Also shown are the models (total, cluster, and combined Galactic fore- and cosmic X-ray background) from the simultaneous fit as 
well as the residuals in terms of standard deviations. There are no strong systematic deviations for any of the regions or detectors. 
Right panel: Same as left but only for bin 3 (7.5'— 1 1 .5'). In this plot the Galactic fore- and cosmic X-ray background models 
are shown separately, and the subtracted particle background is also given (by the thicker data points at the bottom). The cluster 
emission dominates over the other components only in the range ~1— 1 .5 keV and even there only marginally (see text). 



(marginally) around 1 keV until the cosmic X-ray background 
becomes stronger around 1 .5 keV. The particle background starts 
to dominate over the Galactic foreground around 1 .4 keV, over 
the cluster component around 3 keV, and over the cosmic X- 
ray background between 4 and 5 keV. This plot illustrates that a 
cluster temperature determination in this region requires long ex- 
posures and is getting close to the limit of Suzaku's capabilities. 
In case of Chandra or XMM-Newton, the particle background 
would start to dominate over each of the other components at 



much lower energies. The same plot for bin 4 looks quite simi- 
lar, there are just fewer data points and no cluster emission. 

Figure Q] (right) shows projected gas temperature profiles 
measured with Chandra, XMM-Newton, and Suzaku. In the in- 
nermost region (<0.4') the XMM-Newton temperature is con- 
sistent with the range of best fit temperatures obtained with 
Chandra (3.5-7.6 keV). In the region ~l-2', however, Chandra 
gives systematically higher temperatures than XMM-Newton. 
The primary reasons for this are, first, that the XMM-Newton 
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Table 1. Best fit Suzaku cluster temperatures, metallicities, 
and the corresponding statistical uncertainties (90% confidence 
level) of the four radial bins. 



Bin 


Radius/ ['] 


Radius/ [R 2m ] 


kT x /[keV] A /[solar] 


1 

2 
3 
4 


0-3.5 
3.5-7.5 
7.5-11.5 

12-16 


0-0.30 
0.30-0.64 
0.64-0.98 
1.02-1.36 


6 65 +008 

-0.08 

7 C7+O.39 

-0.38 
4 40+0.79 


43 +0 02 

-0.02 

37 +0 07 

"••"-0.07 

25 +0 25 



profile is not corrected for PSF smearing, so the temperatures 
in this region are slightly underestimated; and, second, that the 
Chandra effective area calibration is, apparently, inaccurate, re- 
sultin g in overestimates of the temperature at high temperatures 
(e.g., Sn owden et aHl2008l: ISun et al.ll2008l and L. David's pre- 
sentatior0). The central Chand ra temperatures deter mined here 
are consistent with those from ISanders et all d2005l) . In the re- 
gion ~3-6' Chandra and XMM-Newton give consistent results. 

The Suzaku best fit temperatures and uncertainties are given 
in Tab. [T] The innermost Suzaku bin covers 0'-3.5' and, obvi- 
ously, contains emission from quite a range of temperatures. The 
result of the single temperature fit (6.65+° ^ keV) lies roughly in 
the middle of this temperature range, as expected. The Suzaku 
best fit temperature in the second bin (3.5'-7.5', 7.57^ ^ keV) 
is slightly higher compared to the Chandra and XMM-Newton 
temperatures. The primary reason for this is most likely contam- 
ination from emission from the very bright region <3.5' due to 
Suzaku's PSF (contamination is ~50%, see discussion below). 
Even though A2204 appears quite regular on large scales (e.g., 
ISchuecker et al.ll2001l) . possible deviations from exact spherical 
symmetry of the temperature structure may also introduce some 
scatter in the comparisons. 

For the purpose of this paper, bin 3 covers the most important 
region. It extends from 7.5' to 1 1.5'; i.e., from 64% to 98% of the 
estimated rjoo. The Suzaku best fit temperature (4.49^'gg keV) 
is significantly lower than the temperatures measured further in. 
For Chandra and XMM-Newton, temperature measurements in 
the low surface brightness cluster outskirts are strongly affected 
by the uncertainty in the particle background. We tested the in- 
fluence of the particle background on the Suzaku results by arti- 
ficially changing the nor malization of the N ight Earth spectra by 
a very conservative (e.g jTawa et al.ll2008l) +10% and repeating 
the fitting procedure. We found that the new best fit temperatures 
are consistent with the statistical uncertainty of the standard fit 
(Fig. Hi. This demonstrates that the results obtained here are not 
affected significantly by uncertainties in the particle background 
subtraction. 

We also attempted a temperature measurement with Chandra 
in the region 6.9-9.9'. The Chandra best fit temperature is con- 
sistent with the allowed Suzaku temperature range but the statis- 
tical uncertainty (excluding systematic effects due to the particle 
background) is too large for a meaningful constraint of the tem- 
perature with Chandra. 

4. Discussion 

We showed clear evidence that the temperature of this relaxed 
cluster declines significantly when going from the inner to the 

2 http://cxc.harvard.edu/ccw/proceedings/07_proc/presentations/davic 



outer regions. This is expected theoretically (e.g., iFrenk et all 
1999). To compare, in detail, the result obtained here to predic- 
tions, we overplotted in Fig.[3](left) the average temperature pro- 
file in cluster outskirts as determi ned with hydrodynam ical sim- 
ulations of massive clusters by iRoncarelli et alj d2006l Sample 
A in their Tab. 2). These authors specifically studied the regions 
around the virial radius. Note that magnetic fields and cosmic 
rays were not included in their simulations. For the comparison, 
we used r 200 = 11.75' and kT x (03r 2 oo) = 7.42 keV. The lat- 
ter was calculated by taking the average of the best fit Chandra 
and XMM-Newton temperatures of the bins that include 0.3r2oo 
(shown in Fig. [3] left). We excluded the Suzaku measurement 
because its broad PSF makes an accurate determination at this 
radius difficult. 

The temperature of the inner Suzaku bin shown is slightly 
higher but consistent with the prediction (keep in mind that, 
in general, regions left of bin centers carry more weight in the 
temperature determination than regions right of bin centers, be- 
cause the surface brightness decreases rapidly with radius). As 
mentioned above and discussed in more detail below, this bin is 
strongly affected by PSF smearing, so the actual projected tem- 
perature may be lower. Due to the steep surface brightness pro- 
file, deprojection would likely result in only a minor increase of 
the temperature estimate. The important outer bin is much less 
affected by PSF and projection effects and is slightly lower but 
consistent with the prediction. In conclusion, the observed outer 
temperature profile is slightly steeper but consistent with a drop 
of a factor of 0.6 from 0.3 r2oo to r2oo, as predicted by simula- 
tions. 

We tested illustratively what the improvement in the uncer- 
tainty of a total mass estimate is due to the improvement of the 
uncertainty of the temperature profile provided by the Suzaku 
data. We fitted powerlaws to the outermost XMM-Newton data 
point and the Chandra and Suzaku outermost data points, re- 
spectively, taking into account the upper and lower statistical 
(90%) errors of the latter two. Then we determined a fiducial 
total mass using the best fit single fi model para meters for the 
density profile (from Reiprich & Bohringer 2002) and the pow- 
erlaw model parameters for the temperature profiles, under the 
assumption of hydrostatic equilibrium. Using the Chandra tem- 
perature we found M(< 1840 kpc) = 5.34+J-^ x 10 14 M G . Using 
Suzaku we found M(< 1840 kpc) = 6.18+^jj x 1O 14 M ; i.e., the 
uncertainty due to the temperature profile decreased from 40- 
60% (Chandra) to 10% {Suzaku). So, this represents a signif- 
icant improvement of a factor of ~5 (factor of ~3 when taking 
into account the different exposure times). Note that a full analy- 
sis of the total mass profile will be performed in a more in-depth 
work. 

In the following, we discuss and quantify several system- 
atic effects. We discuss these issues here at some length because 
this is one of the first papers determining cluster properties out 
to ^200, where the surface brightness is very low and system- 
atic effects potentially quite important. In summary, we found 
that the combination of all quantified systematic effects is of the 
same order as the statistical uncertainty in the most relevant bin 
(7.5'— 1 1 .5'). The individual results of the tests are summarized 
in Tab. |2] 

We start with the particle background. The reproducibility 
of this ba ckground compon ent for Suzaku is better than about 
3% (e.g., Taw a et all 120081) . Here, we conservatively assumed 
5% and found best fit temperatures and errors for the third bin 
as 4.66+J}^ keV and 4.32+J} 77 keV. So, for Suzaku, this effect is 
very smail. 
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Fig. 3. Left pane l : Obse rved outer temperature profile compared to profile and scatter predicted by hydrodynamical simulations of 
Roncarell i et al.l (|2006, solid lines). Symbols have the same meaning as in Fig. [TJ For clarity, only the two Chandra and XMM- 
Newton data points are shown that were used to determine &rx(0.3r2oo). The other Chandra and XMM-Newton data points further 
out are consistent with the simulation results. The dashed error bars of the outer Suzaku bin indicate the combined statistical plus 
systematic error range, calculated as described in the text and using the uncertainties given in Tab. [2] Right panel: Contamination 
fractions of the XIS 1 bin 7.5'— 1 1 .5' determined from ray tracing simulations. 75% of the photons detected in this bin also originated 
in this region. The energy dependence is negligible (the lines for all five energies overlap) and the results are very similar for the 
other detectors. 



Next, we discuss the influence of Suzaku's PSF. We 
performed ray tra cing simulations using the xissim tool 
( IIshisakietall 20Q7) using 10 7 photons per energy and detector. 
We assumed Abell 2204's surface brightness profile to follow the 
best double (3 model fit to the Chandra dat a. We followed a pro- 
cedure very similar to the one described bv lSato et al.l (120071) and 
found that the second bin (3.5'-7.5') is significantly (49%) con- 
taminated by emission originating in the cluster center (<3.5')- 
The most important region, the third bin (7.5'— 1 1 .5') is contam- 
inated by 25%; i.e., a relatively small fraction (see Fig. [3] right, 
for the contamination fractions of this bin determined for XIS 1 , 
as a representative example). The energy and detector depen- 
dence of the PSF were found to be small enough to be negligible 
here. The same is true for the fine details of the model for the 
surface brightness profile, since we found very similar fractions 
when repeating this analysis using the best fit double /3 model as 
obtained from ROSAT PSPC pointed data. 



Naively, one could assume that a simple PSF contamination 
correction could be performed, starting with the assumption that 
49% of the emission in the 3.5'-7.5' bin comes from a plasma 
at kTx = 6.65 keV (Tab. [TJ. However, this would not be accu- 
rate, as the Chandra and XMM-Newton data reveal a wide range 
of temperatures at different distances from the cluster center 
(Fig- IB- A detailed correction for the PSF effects, therefore, re- 
quires a large number of ray tracing simulations to be performed, 
using fine Chandra temperature bins, to determine precisely how 
much emission from plasma at what temperature contaminates 
each of the Suzaku bins. A detailed correction for the PSF ef- 
fects is beyond the scope of this paper. Such a correction will be 
performed in a more in-dep th analysis, using th e longer Chandra 
observation of this cluster (Sanders et al. 2009) and the recently 
updated effective area calibration. Here, we follow a rather con- 
servative route and quantify the changes in temperature and un- 
certainty by assuming the temperature for the emission, contam- 
inating the important third bin, to lie in the range 6.0-8.0 keV. 
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Table 2. Systematic uncertainties of the cluster temperature 
measurement in the bin 7.5'— 1 1 .5'. 



Test 


A^r x /[keV] 


1 Particle background (-5%) 


+0.18 


2 Particle background (+5%) 


-0.17 


3 PSF (25% contamination) 


-0.68 


4 Foreground modeling (double thermal) 


+0.61 


5 ROSAT foreground (thermal) 


+0.32 


6 ROSAT background (powerlaw) 


+0.46 


7 ARF from ROSAT surface brightness 


+0.01 


8 Extrapolated cluster emission beyond ?"2oo 


+0.24 



Furthermore, we freeze the metallicities of the third bin to 0.2 
and that of the contaminating emission to 0.3. This results in 
the following best fit temperatures and uncertainties: &Tx(7.5'- 
11.5') = 4.13!°$ keV (for W contami = 6 keV) and kT x (1.5'- 
11.5') = 3.81+°^' keV (for kT cont . dmi = 8 keV). Changing the 
metallicities of the third bin to 0.15 and 0.25 has a quite negligi- 
ble influence on the best fit temperatures. The same is true when 
changing the metallicity from the contaminating emission from 
the second bin to 0.25 and 0.35. In Tab. [2] and the calculation 
of the total systematic uncertainty we use the more conservative 
(£r contami = 8 keV) result. 

We specifically designed the Suzaku observation in such a 
way as to have a cluster free region available to help constrain 
Galactic fore- and cosmic X-ray background directly and locally, 
using the same observation. This ensures that the results are not 
affected by systematic calibration differences between different 
satellites and varying point source subtraction fractions. Still, we 
checked whether the resulting parameter values are in a reason- 
able range. We found that the normalization of the powerlaw 
component expected from deep CXB studies using several dif- 
ferent satellites is lower by factors 0.77-0.62 compared to our 
best fit normali zation. We then fro ze the powerl aw normalization 
to the highest ( iVecchi et al.| [T999) and lowest dRevnivtsev et alj 
120031) normalization we found in the literature, resulting in tem- 
peratures 6.08!g g2 and 6.84^^2 keV, respectively, in the third 
bin. This procedure results in worse fits but note that, in the lat- 
ter case, the temperature is significantly hotter than when leav- 
ing the powerlaw normalization free. Since intensity variations 
of the CXB may be correlated with large scale structure, we 
decided to use the independent ROSAT PSPC observation of 
Abell 2204 and fitted the same fore-/background model in the 
energy band 0.3-2.4 keV to the outskirts. We found that the nor- 
malization is fully consistent (factor 0.81-1.04) with the higher 
powerlaw normalization from the Suzaku data of Abell 2204. 
Moreover, we also found the temperature and normalization of 
the thermal component to be fully consistent with the Suzaku re- 
sults. In particular, for the Galactic emission, we found tempera- 
tures of 0.25+°^ keV (Suzaku) and 0.26+° °| keV (ROSAT), quite 
typical of this foreground component; i.e., there was no indica- 
tion of significant hot cluster or accretion shock emission "con- 
taminating" annulus 4. We, therefore, conclude that our Galactic 
fore- and cosmic X-ray background modeling is adequate. 

Recall that the purely statistical uncertainties on these fore- 
and background components are already included in our quoted 
statistical errors on the cluster temperatures because in our stan- 



dard analysis we performed a simultaneous fit of all the astro- 
physical fore- and background components together with the 
cluster component; i.e, the relevant parameters for these models 
were all free to vary. For the calculation of the total systematic 
uncertainty, we additionally included effects, which relate to the 
modeling of the Galactic thermal emission as well as the spatial 
variation of the fore- and background components, as described 
in the following. 

Several authors used more than one spectr al component to 
model the Galactic foreground emission (e.g JSato et alll2007l: 
ISnowden et all 120081; Benlev & Sheltonl 120081) . Therefore, we 
tried a two component foreground model (phabs*apec+apec) in- 
stead of a single component model (apec) for the Suzaku obser- 
vation. As best fit temperatures for the Galactic components we 
found 0.28 keV for the absorbed component and 0.13 keV for 
the unabs orbed one, in good agr eement with the typical values 
quoted by Snowde n et alj d2008l) . Furthermore, the new power- 
law model normalization (representing the unresolved cosmic 
X-ray background) is only 2.5% lower than in the original fit. 
This change is much smaller than the statistical error on this nor- 
malization. 

We found that the resulting new best fit cluster temperature 
in the most important and most affected third bin (7.5'— 1 1 .5') 
gets a bit higher but not significantly so (5.10 keV). On the one 
hand, this lends confidence to our approach; on the other hand, 
this result is also not too surprising since in our spectral analysis 
we ignored all photons with energies E<0.7 keV, so the cluster 
temperature measurements are by construction less sensitive to 
uncertainties in the soft fore- and background. 

Also, x^. does not change by adding this additional model 
component, showing that no significant improvement can be 
achieved by adding this second thermal component for the data 
under consideration here. 

Moreover, the relative statistical uncertainties stay very simi- 
lar, e.g., for the best fit cluster temperature in the third bin the er- 
rors at the 90% confidence level are-13%+18% (single thermal 
Galactic foreground model) and -15% +21% (double thermal 
Galactic foreground model); therefore, the simple single thermal 
model does not result in a significant underestimate of statistical 
errors. In our systematic error analysis we do include the model 
dependence in the final error. 

Last not least, we tested the double thermal model in the 
ROSAT analysis of the foreground emission. We found very sim- 
ilar best fit temperatures for the Galactic components (absorbed 
0.25 keV, unabsorbed 0.13 keV), although ROSATs poor energy 
resolution results in significantly enlarged errors for the tempera- 
tures of the two thermal components (the two temperature ranges 
overlap; degeneracies cannot be broken, resulting in an unsta- 
ble error analysis). The change in the powerlaw normalization is 
again well within the statistical uncertainty. 

Next, we tested the influence of a possible spatial variation 
of the fore- and background estimates. The regions selected for 
the fore- and background analysis in the ROSAT observation dif- 
fer from those in the Suzaku observation. Therefore, we refit the 
Suzaku data but this time not letting the powerlaw and thermal 
component vary freely but, separately, freezing them to the val- 
ues determined from the ROSAT observation (after correction for 
the different covered areas). The best fit cluster temperatures for 
the third bin do not change significantly. They are 4.81^}^ keV 
(foreground apec model frozen) and 4.95!g ^ keV (background 
powerlaw model frozen). Both effects are included in the total 
systematic error analysis. 

The Monte Carlo ARF calculation using xissimarfgen re- 
quires a priory knowledge of the cluster surface brightness dis- 
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tribution, which we provided using a double (3 model fit to the 
surface brightness measured with Chandra. We tested the influ- 
ence of deviations from this assumed distribution on the clus- 
ter temperature measurements by instead generating ARFs using 
the best fit double j3 model from the ROSAT observation. The re- 
sulting best fit temperature in the third bin was almost unchanged 
(4.50 keV). 

We assumed that the "cluster free" region (the forth bin) con- 
tains negligible cluster emission and argued that the best fit tem- 
perature for the Galactic foreground emission we found supports 
this assumption. Nevertheless, we performed an additional test 
to estimate the contribution of this assumption to the total sys- 
tematic error. To this end, we determined a rough surface bright- 
ness profile using bins 1 through 3 and extrapolated it by conser- 
vatively assuming the surface brightness to drop by a factor of 5 
from bin 3 to bin 4. Furthermore, we assumed a temperature of 2 
keV and a metallicity of 0.2 solar in bin 4. We then included this 
cluster component in the model for bin 4, in addition to the fore- 
and background components, and redid the full simultaneous fit. 
The resulting new best fit cluster temperature in bin 3 is 4.73 
keV; i.e., slightly higher than the original fit result but not sig- 
nificantly so. This rather small change is expected because the 
extrapolated surface brightness of this fourth bin is much lower 
than the surface brightness of the powerlaw background compo- 
nent. Therefore, the new normalization of the powerlaw back- 
ground component is also only slightly lower (7%) compared 
to the original one. A detailed determination and discussion of 
the surface brightness profile and especially its error will be pre- 
sented in a subsequent paper. We do add the effect of the possible 
cluster emission in bin 4 to the systematic error budget. 

Some of the systematic effects described above raise the tem- 
perature (e.g., using a double thermal model for the Galactic 
foreground emission) others lower it (e.g., roughly accounting 
for PSF effects). To estimate the total systematic error as well as 
the combined statistical plus systematic error, we used the fol- 
lowing scheme. We separately added all positive and negative 
errors from Tab. [2] in quadrature, resulting in a total systematic 
error keV. The total systematic error is, therefore, slightly 
larger but of the same order as the statistical uncertainty C^lg 
keV). Then, we added the systematic and statistical errors in 
quadrature, resulting in a combined statistical plus systematic 
error keV (shown as dashed error bars in Fig. (3). 

For the comparison to simulated temperature profiles, the 
uncertainty of r2oo as estimated from the extrapolated XMM- 
Newton density and temperature profiles may be important. For 
instance, using a c luster mean temperature and the relation of 
Evrard et al J (119 96) would result in a larger value for r2oo- Also, 
Arnaud et all (120051) determined r2oo = 2075 kpc using the same 



XMM-Newton data, probably due to the flatter temperature pro- 
file they found. This would make the observed temperature drop 
even steeper in terms of r2oo, potentially resulting in tension 
between observation and simulations. We will redetermine r2oo 
from the cluster mass profile, taking advantage of the informa- 
tion on the outer temperature from Suzaku in the more in-depth 
analysis of this cluster. It is reassuring to note that further corrob- 
oration for the r2oo estimate used here com es from the indepen- 
dent weak lensing analysis of this cluster (IClowe & Schneider! 
|2002|) . yielding r2oo = 1 1.8', in perfect agreement with our esti- 
mate. 

The confirmation of the predictions from hydrodynamical 
simulations for the gas physics in cluster outskirts indicated here 
rests on the analysis of a single cluster. What is required for 
a general confirmation is the analysis of a statistical sample of 



clusters. In the future, this can be performed with dedicated clus- 
ter observations accumulating rapidly in the Suzaku archive. 
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